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Background Of The Invention 

[001] 1. Related Applications. 

[002] This application claims the benefit of U.S. Provisional Application No. , 

filed September 2, 2003, titled "A GPS System." 

[003] 2. Field of the Invention. 

[004] This invention relates generally to signal processing in a radio receiver and more 
particularly to removal of undesired signals received at the radio receiver. 

[005] 3. Related Art. 

[006] Numerous communication systems rely on receiving one or more radio frequency 
(RF) signals. As the available frequency bands become more congested, the numerous radio 
frequencies used by the different communication systems start to interfere with each other. One 
type of interference encountered in the different communication systems is carrier wave 
interference or more commonly called carrier wave (CW) jamming. 

[007] CW jamming is a source of interference in spread spectrum systems, such as CDMA 
cellular telephone systems and satellite positioning systems. Spread spectrum communication 
systems use lower power signals spread across the frequency spectrum and are subject to 
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interference from carrier waves used in other communication systems. The problem of CW 

jamming is further complicated because of the geographical area covered by spread spectrum 

system may include the whole Earth, for example the United States' Global Position System 

(GPS). 

[008] Attempts to eliminate CW jamming signals have occurred by regulating frequency 
use and by adding dedicated circuitry to radio receivers. The added circuitry being included in 
radio receivers often result in additional power consumption and expense. Some examples of 
jamming signal detection approaches and signal performance improvement approaches may be 
found in United States Patents 6,448,925, 6,577,271, 6,606,349, and patent application 
09/562,749. 

[009] Therefore, there is a need for methods and systems for identifying and removing CW 
jamming signals that overcomes the disadvantages set forth above, and others previously 
experienced. 

Summary 

[010] Systems consistent with the present invention provide a receiver that is capable of 
receiving a spread spectrum signal that contains a CW jamming signal along with a weak signal. 
The signal is processed with a crosscorrelator that enables the CW jamming signal to be 
identified, tracked, and reproduced. The replicated CW jamming signal is subtracted from the 
received signal after demodulation, thus enabling the weak signal to be processed. 
[011] Other systems, methods, features and advantages of the invention will be or will 
become apparent to one with skill in the art upon examination of the following figures and 
detailed description. It is intended that all such additional systems, methods, features and 
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advantages be included within this description, be within the scope of the invention, and be 

protected by the accompanying claims. 



Brief Description of the Figures 

[012] The components in the figures are not necessarily to scale, emphasis instead being 
placed upon illustrating the principles of the invention. In the figures, like reference numerals 
designate corresponding parts throughout the different views. 

[013] FIG. 1 is a flow diagram depicting strong signal cancellation in a weak spread 
spectrum signal using crosscorrelation. 

[014] FIG. 2 illustrates a signal-processing diagram for identifying and removing CW 

jamming signals where the cancelled signal of FIG. 1 may be the CW jamming signal. 

[015] FIG. 3 is a block diagram of electrical components of FIG. 2. 

[016] FIG. 4 is a flow diagram of the identifying and removing CW jamming signals. 

Detailed Description 
[017] Turning first to figure 1, a flow diagram 100 depicting strong signal cancellation in a 
weak spread spectrum signal using crosscorrelation is shown. A strong/weak or near/far signal 
isolation provided by a spread spectrum, pseudo random number (PRN) code family such as 
used in CDMA spread spectrum systems is dependent upon the crosscorrelation between the 
various code members of the family. In the case of a satellite positioning system, such as GPS, 
the isolation of two signals at the same frequency (or multiples of the code repetition rate, in this 
case 1 KHz) is about 21 to 23 dB. If the relative strengths of two signals differ by more than 
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this limit, the weaker signal cannot be discriminated using only the spreading code. A method of 

removing the effects of the stronger signal may be applied if the weaker signal is to be tracked. 

[018] The crosscorrelation effect is at its maximum when the relative Doppler frequency 

offset between the relatively strong and weak signals is an integer multiple of 1 KHz in the case 

of coarse acquisition code (C/A) in the GPS signals. A general solution to the problem of 

tracking a weak signal spread spectrum signal in the presence of a stronger spread spectrum 

signal is based on the premise that all aspects of the strong signal's interference can either be 

measured or calculated in order to remove it from the weaker signal. The solution can be 

implemented in any multi-channel receiver having the ability to control a channel's frequency 

and phase as well as selecting the desired spreading code and setting that code's phase position. 

The receiver typically employs two channels, one to track the weak signal and one to track the 

interfering strong signal. However, the channel that is used to track the strong signal is not 

required if the characteristics such as power, code phase and frequency of the strong signal can 

be obtained or accurately estimated by alternate means. 

[019] As shown generally in FIG. 1, the procedure starts 102 with the strong signal being 
acquired 104, such as by tracking the strong signal in a first channel of the receiver. The channel 
provides a measurement of the signal strength of the strong signal along with the phase of the 
carrier signal and the spreading code. Additional channels may be used to track additional 
strong signals (not shown in FIG. 1). 

[020] The code phase of the spreading code of the weaker signal, along with its received 
frequency and signal phase, are predicted 106 based on the 50Hz navigation data code data (D) 
by methods known in the art. A second channel in the receiver is dedicated to receiving the 
compound carrier signal and tracking 108 the predicted weak signal component. 
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[021] The second receiver channel correlates the incoming signal with the second code at 

the predicted frequency and signal phase. The resulting in-phase and quadrature (I,Q) 

measurements contain both the weak signal and the strong signal, each spread by their unique 

code. Correlation by multiplication of the replica code for the second signal, Code2R, with the 

incoming signal yields the product Code2R*(weak2*Code2+StrongX*CodeX+ . . .) where 

weak2 is the power of weak signal 2, Code2 is the actual code for satellite 2 broadcasting the 

weak signal 2, StrongX (X=l, 3, 4, . . . ) is the power of strong signal X, and CodeX is the actual 

code for satellite X contained in the signal. The product Code2R*Code2 is the autocorrelation 

of the received code 2 and the replica code 2. The autocorrelation function has a value of 1 if 

the replica code is aligned with the received code. This crosscorrelation of replica code 2 with 

Strong X (Code2R*CodeX) is next computed 1 10 to be removed from the compound signal. 

[022] Codel and Code2 are both members of a PRN code family, and their autocorrelation 

and crosscorrelation properties are known. It is therefore possible to calculate the 

crosscorrelation of the two codes at their respective phases by simply multiplying each bit of 

Codel by the corresponding (in time) bit of Code2 to produce their crosscorrelation value. 

Since there may be a relative Doppler frequency offset between the two codes, the phase of the 

codes will process past one another over time and create a new crosscorrelation function. For 

the GPS system the greatest delta code Doppler typically encountered is about plus or minus 

9KHz which is equivalent to six code chips per second (1540 carrier cycles per code chip), and 

thus the maximum recalculation rate of the crosscorrelation value is roughly 6 times per second. 

[023] The maximum crosscorrelation occurs at a frequency offset of zero with peaks 

occurring at intervals of 1000Hz. There is an attenuation of the crosscorrelation as the 

frequency offset moves away from zero. This attenuation follows the well known sin(x)/x 
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curve. If 10ms measurements are used for tracking or acquisition, the attenuation factor would 

be equal to sin(Afreq*7t/100 Hz)/(Afreq*7c/100 Hz). This produces an attenuation of -lOdB at 

about a 75Hz delta frequency. Other local peaks in the sin(x)/x curve (i.e. locally minimum 

attenuation) occur at 150Hz and 250Hz with attenuations of -13.5dB and -18dB, respectively. 

This implies that for a desired strong signal suppression of lOdB, only the first lobe of the 

sin(x)/x function need be considered; however, should additional suppression be desired, the 

entire curve may be considered. 

[024] The next step entails computing 112, for each strong signal, the product of the strong 
signal amplitude and the calculated frequency and time domain (code phase) crosscorrelation. 
The weak signal is finally extracted by subtracting 1 14 this product from the compound signal 
and processing is complete 116. The weak signal thus extracted is subsequently processed in the 
receiver circuitry as known in the art. 

[025] The in-phase (I) and quadrature amplitude (Q) of each strong signal is obtained by 
measurement in each strong signal's own individual receiver channel or by estimation through 
independent means. Because the strong signal is being actively tracked by the receiver's phase 
lock loops, the phase of the strong signal is presumed to be near zero radians and thus nearly all 
the signal power is in the in-phase portion. 

[026] A signal comprising a strong signal SI modulated with a first code Codel summed 
with a weak signal w2 modulated with a second code Code2 produces (Sl*codel+w2*Code2). 
The sum of the two signals is correlated with a replica of the second code Code2R to produce 
£{Code2R*(Sl*Codel+w2*Code2)}, where the sum E includes all chips of the PRN code used 
to modulate the weak signal w2. The autocorrelation of a code with itself is 1 so the preceding 
equation can be rewritten as S{Sl*Codel*Code2+w2}. It can see that in order to obtain w2 
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Sl*Codel*Code2 is removed. Since we know Codel and Code2, we can easily calculate their 

crosscorrelation. This leaves us to estimate the value of SI that can be done by independently 

tracking the strong signal on a separate channel, or by any other convenient means. This 

computed value of Sl*Codel *Code2 would be sufficient if the strong signal SI and the weak 

signal w2 where at the same frequency. The two signals are received at different frequencies, 

however, due to the Doppler effect as well as other factors enumerated previously. 

[027] We know that strength of the crosscorrelation varies with the difference between 

these frequencies in a sin(x)/x relationship. We may therefore calculate an attenuation factor 

based on the difference in frequency between the strong and the weak signal and apply it to the 

computed crosscorrelation. Furthermore, if more than one strong signal is present, an 

attenuation factor is computed for each strong signal. 

[028] The code dependent portion of the crosscorrelation factor is computed from the 
known relative states of the PRN code generators to predict the crosscorrelation between a 
strong signal of unit power and zero frequency offset, and a weak signal. This factor is 
multiplied by the amplitude of the corresponding strong signal and adjusted for frequency 
attenuation before it is subtracted from the composite signal. 

[029] The various Gold codes used to modulate the PRN signals are all derived from a two 
code sequence Gl and G2 where the bits of the two code sequences are combined through an 
XOR operation after G2 has been offset some number of bits relative to G 1 depending on the 
Gold code selected. It is known that an XOR operation using binary numbers is mathematically 
equivalent to multiplication of ±1. This allows expressing the equations below in term of 
products of ±1 while in reality the implementation could be with binary numbers with XORs. 
[030] The correlation between two C/A codes can generally be expressed as: 
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[031] ESatlGl(I)*SatlG2(I)*Sat2Gl^ 

[032] where 

[033] I=Summation index ranges from 0 to 1022. 

[034] SatlGl(I)=Value of satellite l's Gl coder chip at state I. Possible values are ±1 . 

[035] Satl G2(I)=Value of satellite 1 f s G2 coder chip at state I. Possible values are ±1 . 

[036] Sat2Gl(I)=Value of satellite 2's Gl coder chip at state I. Possible values are ±1 . 

[037] Sat2G2(I)=Value of satellite 2's G2 coder chip at state I. Possible values are ±1 . 

[038] offset=time difference between the satellite 1 and 2 in units of chips A0*Phase 

change per chip between satellite 1 and 2 in radians. 

[039] It should be noted that when the difference I-offset is less than 0, 1023 is added to the 
difference to maintain the value in the range of 0 to 1022. In other words, the domain of the 
functions returning coder chip states is limited to the range of 0 to 1022. 

[040] The computation time required to compute the 1023 bit-by-bit correlations can be 
accelerated by making use of standard controller instructions that perform 8, 16 or 32 bit-wise 
XORs with a single controller instruction. The following will demonstrate the method of 
computing eight chips in parallel. Those skilled in the art will immediately recognize that the 
scheme can be easily modified to accommodate some other convenient number of bits per 
controller XOR operation. Such an approach may employ a single 1024 bit correlator match 
filter that performs the correlation for 1023 G1/G2 states concurrently. 

[041] The 1023 states of Gl and G2 may be stored linearly in permanent memory. Thus, it 
is possible to quickly gather 8, 16, 32 or some other convenient number of bits with a single 
controller load instruction by computing the address of the desired chip and the shift required to 
align it. Thirty-two bits is a particularly convenient number because 31 divides 1023 evenly. If 
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a single 1024 bit correlator match filter is employed, the 32 bits of the PN code are derived at a 

time from the G1/G2 states and loaded into a 1024 wide reference shift register. The current 

implementation thus reads 32 bits at a time and uses 3 1 of them at a time for each of 33 intervals 

that span the 1023 chips of the C/A code. The 31 bit sums are broken into four parts of 8, 8, 8, 

and 7 bits, and each 7 or 8 bit sum is multiplied by e" jAGI where I changes by 7.75 chips for each 

part. The form of the sum is: 

[042] E(d A0I *31*I(SatlGl(I^ 

J-offset))+e- jAe(,+3K7 75) *S(Satl Gl (1*3 RJ+8)*SatlG2(I*31+J+8)*Sat2Gl(I*3HJ+8-offset) 
* Sat2G2(I * 3 RJ+8-offset))+e' jA0(I * 3 1 + 1 55) *S(SatlGl(I*31+J+16)*SatlG2(I*3HJ+16) 
*Sat2Gl(I*31+J+16-offset)*Sat2G2(P 
*Sat 1 G2(l*3 1 +J+24)*Sat2Gl (1*3 
[043] where 

[044] I=Outer index ranges from 0 to 32 

[045] J=Inner index ranges from 0 to 7 for the first three sums and from 0 to 6 for the 

last sum. The inner sums are computed in parallel by using a 32 bit word that 
contains all 31 bits and using bitwise XOR to perform the multiplications and 
shifting and adding to sum the 1 bit products. 
[046] Note that all of the multiplications of the Gl and G2 codes in the above equation are 
implemented by bit-wise XOR instructions. The above algorithm is in error by at most -17 dB 
from an exact computation, and requires about 6000 controller operations to complete. 
[047] Periodically, as needed, the code dependent crosscorrelation factors are computed for 
all strong and weak signal pairs with small frequency differences, i.e. frequency differences that 
could cause strong-weak crosscorrelation interference. In the current implementation strong 
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signals are those with C/No >40 dB and weak signals are those with C/No >30 dB. Because 

1 0ms integrations of I and Q measurements are used by the code and phase tracking loops, the 

maximum "significant" frequency difference (modulo 1000Hz) is 90Hz. In the preferred 

embodiment the code dependent cross correlation factor for each possibly interfering pair of 

signals is computed for each of the measurements that might potentially be used by the tracking 

and signal processing algorithms. For example, if early, punctual and late measurements are 

used by the tracking loops, the correlation factors for each of these code alignments is computed 

and stored in the tables. 

[048] These tables only need be updated at a 10Hz rate because the maximum Doppler 
difference is less than 9KHz or less than six chips per second. In addition to maintaining the 
crosscorrelation table, the frequency attenuation of the crosscorrelation due to the frequency 
difference is computed at the 10Hz rate. The attenuation can be expressed as: 
[049] Frequency Attenuation=sin (AF mod 1000*tc/100)/(AF mod 1000*tt/100) 

[050] where 

[051] AF=Frequency difference between a strong and weak signal in Hz 

Mod=modulo offset to give a range of -500 Hz to +500 Hz 
[052] The attenuation only needs to be recomputed if the frequency difference changes by 
more than 5Hz. 

[053] An estimate of the phase and amplitude of the strong signal is required to remove the 
crosscorrelation. The method used in the preferred embodiment is to track the strong signal on 
its own dedicated channel and collect the I, Q measurements output over the exact same interval 
that the weak signal I, Q samples are taken. The known phase and frequency of the replica 
signal that is used to track the strong signal is an excellent approximation of the actual phase and 
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frequency of the strong signal. Furthermore, because the strong signal is in phase lock, the 

magnitude of the I measurement provides a good approximation of the amplitude of the strong 

signal. Finally, the bi-phase modulation of the strong signal data bits D may cause the phase of 

the strong signal to rotate 1 80 degrees whenever the data bits transition from a 1 to 0 or from a 0 

to 1. In the current implementation, the phase of the strong signal is corrected by adding 180 

degrees to the phase of the replica signal whenever the sign of the I measurement for the strong 

signal is negative. 

[054] Every 10ms a new set of I, Q correlation data is available from the channel assigned 
to track the weak signal. The tables of crosscorrelation factors are checked to predict the 
presence of any interfering strong signals. If strong signals are predicted, the following 
subtraction is performed to remove the strong signal crosscorrelations: 
[055] FirstCodeOffset=WeakCodeState-StrongCodeState-StrongDoppler*AT- 

TableEntryOCodeState 
[056] DeltaPhase=WeakCarrierPhase-StrongCarrierPhase-StrongDoppler*AT 

+DeltaKHz* StrongCodeState 
[057] FirstPhase=FirstCorrelationPhase+DeltaPhase 
[058] SecondPhase=SecondCorrelationPhase+DeltaPhase 
[059] FirstMag=FirstCorrelationMag+FirstCodeOffsetFraction*StrongI 

*FrequencyAttentuation 
[060] SecondMag=SecondCorrelationMag*(l-FirstCodeOffsetFraction) 

* StrongI * Freq uency Attenuation 
[061] CorrectedWeakIQ=WeakIQ-FirstMag*e jFiretPhase -SecondMag*e 
[062] where 
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[077] 



WeakIQ=IQ correlation from the weak signal's channel. 



[078] CorrectedWeakIQ=IQ correlation corrected for crosscorrelation from the strong 



signal. Corrected WeaklQ is computed for the early, on time, and late correlators 
by shifting the FirstCodeOffset appropriately, such as by half a chip each. These 
modified correlations are then used normally in the carrier and code tracking 
software for the weak signal. The algorithm attenuates the crosscorrelation by at 
least 10 dB without attenuating the weak signal, and is repeated for each strong 
signal that may be interfering with the weak signal. 



[079] While FIG. 1 has described crosscorrelation to cancel a strong signal having a PRN 
from a weak signal having different PRN, the procedure may be modified to cancel a CW 
jamming signal and enhance a weak signal having a PRN. In FIG. 2, a signal-processing 
diagram for identifying and removing CW jamming signals where the cancelled signal of FIG. 1 
may be a CW jamming signal is shown. The receiver 200 receives a spread spectrum signals at 
antenna 202. The crosscorrelator 204 is placed in a mode to identify CW jamming signals. In 
that mode, the crosscorrelator employs all ones for a PRN code. The crosscorrelator 204 may 
have a signal processor 206 and a match filter 208 that despreads the spread spectrum signal. 
The resultant signal is the CW jamming signal (the strong signal) and is tracked in the tracker 
210. An example of the tracker 210 in the current implementation may be a phase lock loop 
circuit. 

[080] The identified CW jamming signal is then used to generate a replica CW jamming 
signal 212. The replica CW jamming signal may be created using a signal processor 214 and a 
match filter 216. In other implementations, other approaches to generating a desired signal may 



13 



U.S. Express Mail No.: ET972568215US PATENT 
Filing Date: October 20. 2003 Docket No. ST02010USU (246-US-U1) 

[063]. WeakCodeState=Code state of last output to the weak signals channel 



[064] StrongCodeState=Code state of last output to the strong signals channel 

[065] StrongDoppler=Doppler of last output to the strong signals channel 

AT=Difference in time between outputs to the weak and strong channels 
[066] TableEntryOCodeState=€ode state difference of the first element of the 

crosscorrelation table. 

[067] WeakCarrierPhase=Carrier phase angle of last output to the weak signal channel 

[068] StrongCarrierPhase=€arrier phase angle of last output to the strong signal 

channel 

[069] DeltaKHz=Nearest integer multiple of 1 KHz of the difference between the weak 

and strong channels Doppler. In units of KHz. 
[070] FirstCorrelationPhase=Phase entry in the crosscorrelation table for the chip 

indicated by FirstCodeOffset 

[071] SecondCorrelationPhase=Phase entry in the crosscorrelation table for the chip 

indicated by FirstCodeOffset +1 chip. 

[072] FirstCorrelationMag=Magnitude entry in the crosscorrelation table for the chip 

indicated by FirstCodeOffset. 

[073] SecondCorrelationMag=Magnitude entry in the crosscorrelation table for the chip 

indicated by FirstCodeOffset+1 chip. 

[074] FirstCodeOffsetFraction=Fraction of a chip in FirstCodeOffset. 

[075] StrongI=Absolute value of I correlation from the strong channel. 

[076] FrequencyAttenuation=Attenuation due to frequency offset. 
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be employed. Examples of other approaches include, but are not limited to voltage controlled 

oscillators, digital signal processing, analog signal processing. 

[081] The spread spectrum signal received by receiver 200 is then processed by a 
crosscorrelator 214 with dispreading codes rather than the all ones. The received weak signal in 
the spread spectrum signal is demodulated by the signal processor in the crosscorrelator 218. 
The generated CW jamming signal is canceled from the weak signal by a canceller 220. The 
resultant signal is then tracked by a tracker 222. The CW jamming signal that is being tracked 
by tracker 210 is further processed to remove the desired weak signal from the tracked CW 
jamming signal resulting in a clean replica of the unwanted signal on the demodulation process 
of the desired signal in block 224. 

[082] The identification of the CW jamming signal has been described as occurring with 
three crosscorrelators 204, 218 and 212. But in practice, only one crosscorrelator may be used 
allowing the crosscorrlator to distinguish between PRN codes and to identify and remove CW 
jamming signals. Thus, the generation of a replica of the unwanted signal (CW jamming signal) 
occurs during the demodulation of the wanted signal. 

[083] Turning to FIG. 3, a block diagram of electrical components of FIG. 2 is shown. The 
spread spectrum signal containing the weak signal and the CW jamming signal are received at 
the receiver 200 via antenna 202. The spread spectrum signal is demodulated and the CW 
jammer signal is filtered by the jammer filter 302. The resulting CW jamming signal is tracked 
by jammer tracker 304. The tracked jamming signal is then replicated by a jammer replica wave 
generator 306 as a replicated CW jamming signal. The CW signal has the characteristic of a 
constant phase, therefore the replicated CW jamming signal is scaled and rotated 308 to an 
appropriate phase to be subtracted from the demodulated weak signal. The weak signal is 
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demodulated by a demodulator 310. The phase and magnitude replicated CW jamming signal is 

then subtracted from the weak signal from the demodulator 302 by a signal canceller 312. The 

resulting weak signal is then tracked by a track control circuit 314. The track control circuit 314 

outputs the weak signal to the demodulator 310 and two signal combiners 316 and 318. The 

weak signal is subtracted from the CW jammer signal by the signal combiner 316 and the 

resulting signal is used to scale and rotate the replicated CW jamming signal. Similarly, the 

weak signal is subtracted from the CW jamming signal by signal combiner 318 in order to 

provide a more accurate CW jamming signal to be replicated by the jammer replica wave 

generator 306. 

[084] Even though the current implementation has been described with CW jamming 
signals in general, one skilled in the art would recognize that specific known signals such as IS- 
95 pilot signals, Galileo signals, or other known carriers may be identified, tracked and remove 
from a demodulated desired signal. Unlike other known methods of identifying and removing 
CW jamming signals prior to signal processing or at the front end of a receiver, the current 
approach identifies and removes the CW jamming signal at the back end. This is advantageous 
when desired signals are received at low signal strength, such as in CDMA or satellite 
positioning system. The signal strength is not sufficient to support preprocessing of the signal 
while maintaining the desired weak signals. 

[085] Turning to FIG. 4, a flow diagram 400 of the identifying and removing CW jamming 
signals is shown. The procedure starts 402 with a receiver 200 receiving a spread spectrum 
signal 404. The signal may be filtered and the CW jamming signal is identified 406 using a 
crosscorrelator. Once the CW jamming signal is identified 406, it is tracked 408. The tracked 
CW jamming signal is then replicated 410 by a jammer replica wave generator 306. The 
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replicated CW jammer signal is then subtracted from the received signal 412. The CW jamming 

signal has been removed from the received signal and processing is complete 414. 

[086] The foregoing description of an implementation has been presented for purposes of 

illustration and description. It is not exhaustive and does not limit the claimed inventions to the 

precise form disclosed. Modifications and variations are possible in light of the above 

description or may be acquired from practicing the invention. For example, the described 

implementation includes software but the invention may be implemented as a combination of 

hardware and software or in hardware alone. Note also that the implementation may vary 

between systems. The claims and their equivalents define the scope of the invention. 
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